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cross sectional area of inlet tube 

X-area of conical permeameter at distance x from the 
base 

coefficient of Porchhe inter ' s equation, 
constant for conical permeameter 
constant 

1/S, S = Coefficient of micro roughness 
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(for nonlinear flow experiment) 
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acceleration due to gravity 
wall effect 
head loss 

head difference between inlet and exist at any 
instant 

heads of water at inlet of the conical permeameter 
for time interval t^ and tg- 

hydraulic gradient 

hydraulic gradient m a element at distance x from 
base . 

coefficient of permeability 
specific permeability 
height of conical permeameter 
length of pore channel 
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percentage porosity 

hydraustatic pressure 

discharge through conical permeameter 
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constant depending on shape and surface roughness 
of gravel. 


Reynold's number given by 
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= hydraulic mean radius 

= bottom and too radii of conical permeameter 
= taper of the sides of conical permeameter 
= velocity m x and y direction 
= seepage velocity 
= discharge velocity 
= dynamic viscosity of fluid 
= kmametic viscosity of fluid 
= mass density of fluid 
= weight density of fluid 
= constant. 



ABSTRA CT 


Bor flow "through soils, Darcy's law has been 
generally found to be valid. However, It is so only when 
viscous forces predominate over inertial forces. Incase 
of flow through coarser material like sand and gravel the 
inertial forces and viscous forces are of the same cn'deir of 
magnitude and deviations from the Darcy’s law are - observed. 

In such cases the equation suggested by Forchheimeter 
(i = av + b v ) or Missbach (1 = Cv 1 ) are used for relating 
the hydraulic gradient with velocity. In both the equations, 
the coefficients (a, b) and (C, m) depends on the physical 
characteristics of the porous media and the flow conditionsD 


, A review of the Hiterature shows that considerable 
theoretical and experimental work has been done on non- 

r 

linear flow problens.t/{Bu t no efforts seems to have been 
made to relate the coe f Cicienfcs a and b of the Forchheimer 
equation, specially the coefficient ! b' with the physical' 
characteristics of the media. Therefore^Jhe present 
study is mainly intended to relate the coefficient ! b 1 
with size, porosity, shape and roughness of the gravel 
Since past available data are insufficient for this study 
therefore Experiments were carried out on the gravel m the 
range of sizes from 0.313 to 4.62 cm and glass balls 
from 1.3 cm to 2.5 cm. ' 



Results shows tb-’t 'b ! varies mversly with 


size and porosity of gravel. Two relationships have 
been established between the coefficient r b 1 and the 
noomal size and another between coefficient 'b ! and size, 
porosity, surface roughness and shape of the gravel. A 
new non-dimensional p^r^eter has been developed and its 
variation with porosiry has been presented in the graph. 

Also, the variation between Reynold^ number and friction 

/ 

factor has been shown m another graph, y 


J'k conical permeameter has been developed so that it 

s 

generate a range of hydraulic gradients, sands and gravels 

I 

can be tested m thishsucli that an average value of K 
can be obtained for the range of gradient existing. 
Analytical expressions are obtained for head loss and dis- 
charge velocity through the conical permeameter. Typical 
results for K of fine s^nd obtained from the conical and 
cylindrical permeameter are found to compare well^j 



CHAPTER 1 


INTRODUCTION 

Problems concerning flow of water through sous 
have become one of the most interesting and important 
field of study among the Civil Engineering subjects. This 

is because of new achievements m the fields of Soil 

any 

Mechanics and Hydraulics Engineering. Whenever huge 

structures like dams and other hydraulic struct ares aro 

constructed, one meets with problems like uplift pressure, 

seepage, stability of slopes and many problems related to the 

flow of water m soils. In ground water engineering one 

comes across flow of water m to wflls, trenea?es, etc. 
m 

Much attention is being paid towards such problem. 

In the past as early as 1856 Henry Darcy first 
studied the problem of water flow m compacted sands 
experimentally and put forward the linear relation ship 
between head loss and discharge velocity, which is popularly 
known as 'Darcy’s law. Later on other investigators such 
as Dupuit ,Povlovsky etc. applied the Darcy's law m various 
seepage problems. 

However it has been recognised that there are two 
principal types of flow m soils. 

1. Laminar flow: In which fluid particles travels 

along definite paths which never cross each other. Eor this 


flow Darcy law is valid. 
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2# Turbulent Flow: In which paths are irregular, 

twisting, crossing and re crossing each other at random. 

For the most cases of flow through soils of 
practical interest, the flow velocities and particle size 
are small enough for laminar conditions to exist and 
Darcy's law hold good. To discuss the fundamental approach 
to discribe flow through porous media, it is necessary to 
mention the approach suggested by Philip (1957) and 
developed by Batson (1965). In this, use i^ made of 
Kavier - Stokes' equations of motion which are derived 
from classical hydrodynamic al principles. For newtonian 
fluids, steady motion, incompressible flow and neglecting 
body forces, it can be sho wn (Kay - 1957) that for a two 
dimensional medium. 
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where u,v = velocity m x,y directions 


r mass density 


p = hydrostatic pressure which would occur with 
equilibrium under gravity at the point 
considered. 

v = Kmamatic viscosity 


In these equations the first term represents the 


viscous forces and the right hand side the inertial terms. 



It should be emphasized, however, thar Navier - Stoke 's 
equation using mean flow values are restricted m use of 
laminar flow condition and on the basis of classical 
hydrodynamics, it follows that Darcy's law is at least 
only a good approximation at low Reynoulds number 
since the inertial terms are always present. 

As the flow velocities, particle size and Reynold's 
number increase and the laminar conditions do not exist, 
Darcy’s law becomes increasingly m -accurate and another 
flow law is required. The reason for the non-lmearty in 
the increasing significance of inertia forces as recognised 
by Tek (1957), Green and Duwez (1951), Hubbert (1956), 
Watson (1963). Tek and 0 'Neill (1965) have suggested three 
general regimes of flow through porous media. 

1. Darcy flow regime: In this viscous foi'ces 

predominate, inertial forces are insignificant and Darcy’s 
law is closely follow! . 

2. Laminar non Darcy flow regime * At higher 

Reynold numbers where inertial forces become significant 

the 

as well as viscous forces flow may still be laminar, but 
non- linear and Darcy’s law is no longer obeyed. 

3. Turbulent flow regime: At very high Reynold's 
number, inertial forces become large as compared to viscous 
forces, turbulence sets m and is initially localized m 
large voids, fractures, cavities etc. 
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The range of validity of Darcy’s law has been a 
matter of discussion. Watson (1963) fixes this limit 
between 0.2 to 5.0 while Sohei&igger (I960) fixes as 0.1 
to 75.0. This large variation m upper limit is because 
of the factors like, size, shape, roughness and porosity 

not 

of particle are considered. 

As obvious from category-2~lammar non Darcy 
regime, the hydraulic gradient (l) is no longer directly 
proportional to the flow velocity (v). Therefore some 
relationship between (l) and (v), other than Darcy's law, 
is necessary ior such flow conditions. Forehlieimer (1901) 
and Missbach (1937) have suggested a polynomial and 
exponental forms of relationships between (i) and (v). 

The relationships are : 

2 

i = av + bv (For cheater) (1.3) 

l = cv (Missbach) (1.4) 

where l is hydraalic gradient , 
v is macro valocity , 

a and b are coxistants which depend on size, shape, 
surface roughness of the gravel pack and the flow conditions. 

c and m are also constants depending on flow conditions. 
By analogy with pipe flov/, it could be expected the 
exponent ‘m ' will vary from 1.0 for laminar flow to 2.0 
for fully turbulent flow. 



Nonlinear laminar flow and turbulent flow through 
porous media is encountered m several field situation, e.g. 

a - flow through rock fill dams 
b - flow through filters 
c - flov- is areas adjacent to a well 
d - flow m rock fissures. 

A review of literature (Chapter -2), shows that 
considerable work has been done on nonlinear flow problems. 
The experimental work has been mainly carried out to 
determine the coefficients a, b, and c, ra of Forchhoimer ’s 
Missbnch equations and no attempt has been made to relate 
these coefficient specially a, b with the physical 
characteristics of the porous media. Therefore m this 
work, Chapter-3 is mainly intended to study the coefficients 
a and b of Forchhoimer ’ s equation with a view to relate 
them to the sm: , porosity, roughness of the media. 

In chapter 4 an attempt has been made to develop 
a permeameter to generate a range of hydraulic gradients 
and flow velocities. The properties determined from such 

Of K 

a permeameter could be expected to give average values for 
the range of gradients or velocities generated. Analytical 
expressions useful for this study are given there m. 

Concluding remarks are given m chapter 5 and 
are followed by a set of references. 



CHAPTER 2 


LITERATURE REVIEW 


2. 1 GENERAL : 

A brief review of literature on nonlinear flow 
through, porous media is presented here. The entire work 
can be divided into two categories. Thus, 

(1) Theoretical Work: 

(a) Study of the upper limit of Darcy’s law, 
and determining the laws governing high 
velocity flow through porous media. 

(b) Application of the above theory to the various 
field situations and relating the flow 
characteristics of the porous media to the 
coefficients of the governing laws (Head loss 
equations ) . 

(2) Experimental Work: 

Determination of the values of the coefficients of 
the head loss equations, study of onset of the turbulence, 
determination of the v arious limits of flow regimes, study 
of effect of convergence of stream lines and wall effect, 
study of the relation between friction factor and Reynolds 
number, study of seepage problems by electrical analogue. 
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2 . 2 THEORETICAL WORK; 

(a) Hubbert (1940) reports that Darcy's experiments 

were carried out on a well compacted filter medium con- 
sisting mainly of sand but with proportion of large 
particles. Host of the tests were confined to linear 
flow but he recognised that there was an upper limit to 
the validity of this law. He defined the limit m terms of 
velocity approximately 10 to 11 cms/sec. for this material. 
He has not mentioned the size and energy of compaction 
given for his material. 

Taylor (1948) calculated that for a unit hydraulic 
gradient the maximum diameter of the uniform sand m which 
the laminar flow will occure is 0.5 mm. 

The accepted method of expressing the limitation 
of Darcy's law is to deter fine the Reynold's number at the 
point where departure from the linear relationship is 
first noted. Watson (1955) has carried out a computer 
solution of the Havier Strolce equations for two dimensional 
flow to determine the values of Reynolds number at which 
the permeability becomes noticeably affected by inertial 
terms. The results shows a rapid increase m the inertial 
effect between Reynolds number of 0.2 to 5.0. Scheidegger 
(i960) and Todd (1959) quoxe the ran^e as from 0.1 to 75.0 



and 1,0 to 10 respectively. 

Thus the large variation m the Reynolds number 
reflects the defficiency of the usual form of parameter 
m accounting the variables involved in flow, such as 
size, shape roughness and porosity of the particles. 
Because of this deficiency it was felt necessary to do 
the experimental work m the context of the present study. 

The generally accepted form of head loss equation 
m non Darcy regime, and which has been suggested by 
Porehheimer (1901) is as follows 

i = av + bv 2 (2.1) 

In which a and b are constants depending on the 
physical characteristics of the media and flow conditions 
such as Reynolds number. This equation has been experi- 
mentally investigated by Volker (1968), Ahmed and Sundada 
(1969), O’Neill, Parkin, (1970), Todd and Tyagi (1970), 
Ranganadha Rao and Suresh (1970). By deduction from 
Ravier - Stoke ’s equation Ahmed and Sundada (1969), 
showed that the governing head loss equation is of the 
form suggested by Porchheimer, 

Another form of equation suggested by Missbach 
(1937) is: 


i 


( 2 . 2 ) 



y 


In which c is constant depending on media 
and flow condition and m is an exponent which is 1.0 
for laminar flow and 2.0 for complete turbulent. Above 
equation has Doen supported by Aaandakrishnan and 
Varadara^alu (1963), Parkin, Trollope, lowson (19 66). 


Wilkins (1955) using wide range of aggregate 
sizes, derived the ecuation of the form: 


V v = C 1 * 


(2.3) 


where V ' — seepage velocity 

n 1 = (1/m) 

r* = hydraulic mean radius 

fj. = dynamic viscosity of fluid 

, a^, are constants. 

The value of varied from 32.9 for crushed 
stone to 46.5 for marbles m inch units. 


It may bo noted that Irmay (1958) has remarked 
that although for any type of flow, the nonlinear acceleration 
terms m the Havier - Strokes equation may be neglected at 
small velocities it is only for staight pipes that these 
terms are identically zero at higher velocities. Por 
flow through porous medn their presence m the equation 
will have the effect of causing the resistance - velocity 
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relation to be nonlinear. 

2.2(b) lihmed and Sundada (1969) derived the 

equation from the Navier - Stroke's equation^ of the form 

/X *\ O 

i = — — v + — * — v (2*4) 

/ J gK gfCK 

on comparing this equation with Fo rchhe liner *s eqn. (2.1), 

it has been found that 

, 1 1 

a = rr and b = — , where 

fS K gfCK 

C = constant 

g = acceleration due to gravity 
K = coefficient of permeability 
M = dynamic viscosity 

f = mass density of water. 

Izbash and Lelleeve (1971) derived the formula 
of the form 

l = — 1~— v + ~4 y2 (2.5) 

R^g 05 R 

on comparing with Forchheimer r s equation (2.1) gives 

a = ~ 4 ~ an( 3- b = — \ 

R g Cq R 

where 0 = 1/S , S = coefficient of microroughness. 

R = hydraulic me^n radius 



j3 = coefficient depending on cross sectional shape 
f> = dynamic viscosity 

Izbash and Khaldre (1959) studied the problem 
of turbulent seepage through fill material. Author 
has presented the results in. form of formulae, tables 
and graphs for the seepage through partially or fully 
submerged rimer closures. These closure were made by 
simply dumping tu.e gravels or cribs and which were 
triangular m c ross-sect ion a nd having large voids. 

Me Corquodale and Ng (1969) applied the finite- 
element method to solution of non larcv's flow with a 
free surface problems. 

0 ’Kat (1969), Volkor (1969) also applied the 
finite element method to study nonlinear flow through 
rock fill dam with sheet piles. 

Madhav and Subramanyn (1971) studied the nonlinear 
flow problems arising from confined and unconfined flow 
into wells and trenches. 

Subramanya, Valsangkar and Rao (1971) have 
studied the unconf ined flow over an impervious boundary 
and flow towards a fully penetrating well by using 
Forchheimer 1 s equation. 
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2.3 EXPERIMENTAL WORE: 

Many investigators Rave carried out the 
experimental work and came out with interesting and 
useful results which are presented here m brief. 

Todd (1959) found m his experimental investigation 

that m case of flow through porous media, a particular 
of 

type relation was exhibited between Reynolds number and 

friction factor. He used the formula F., = — to 

1 2 ^ v 2 

calculate friction factor m which d = size of grovel, 

= mass density of water and v is the velocity of 
wat er . 


Parkin (1963) made extensive studies on a, model 
dam m 4' wide flume and found the method of gradually 
varied flow is satisfactory to determine the water surface 
profile . 

Anandakr js hnan end Varada^lu (1963), Dudgeon 
(1964) obtained the values of c and 1 of Missbach equation 
and the values of a and b Forchheimer *s equation for 
different grades of sand and different sizes of gravels 
respectively. 


Dudgeon (1966) obtained the coefficients C and m 


of Missbach eauation. On the basis of his results he 

2 

developed an expression Fgi = F 2 2 “j g » in which Pg 
is friction factor given by ~ , ds is the equivalent 
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spherical dia. , v is velocity, g is acceleration due to 

* 

to gravity and is the length ox channel, 

* 

Dudgeon (1967) determined the wall effect m the 
permeameter. He used Biissbach equation for his investi- 
gations. It was found the exponent m was 1.81 for (wall) 
outer zone and 1.82 for inner zone; the velocity ratio 
between wall zone and inner zone was 1.06 which shows 
that permeability m the wall zone is higher than m the 
inner zone, tho ratio of coeCficient C between inner 
zone and wall zone was approximately 4.0. 

Robmcurtis (1967) has studied the flow through 
rock fill banks ->nd analysed the stability of slope against 
seepage . 

Wright (1969) carried out experiments to detect and 
measure turbulence and to observe the effect of convergence 
of stream lines 021 Reynolds number vs friction factor graph. 
It was found that flow could be classified m foolr regimes. 

1. A laminar regime m which velocity remains unchanged 
at any moment and Darcy's la?/ is valid. 

2. A steady inertial regime m which v elocity remains 
unchanged but Darcy's law is not valid. 

3. A turbulent transition regime in ^hich velocity 
fluctuates at any point with increasing but 
regularly and head loss approximately depends on 
square of velocity. 
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4. A fully turbulent regime m which all parts of flow 

are turbulent , velocity fluctuates randomly about 
a mean, The bead loss now depends on square of 
welocity . 

Volker (1969) obtained the values of a,b and c, m 
of Forchheimer ' s rnd Missbach equations respectively by 
lab studies of model of gravel banks m an open flume. 

Ahmed and Sundada (1969) showed on the basis of the 
results of 18 tests that friction factor (F) and Reynolds 
number (Re) bear a relationship 

F3 = ^ + 1.0 (2.6) 

He used forchheimer 's equation for the investigation. 

Ruble (1969) made extensive study on nonlinear flow 
problems and preseated an electrical analogue model "or the 
solution of such problems. 

Ranganadha R .90 and Suresh (1970) have investigated 
the values of coefficient a and b for Forchheimerfe 
equation for different sizes of gravels. It was concluded 
the relationship gi ( /en by 


1 

Re 


+ 1.0 


(2.7) 
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Proves only the validity of Forchheuner 's equation but 

can not prove the correctness of the values a and b calculated 

from hydraulic measurements. 

Ward (1970) carried out experiments and obtained 
a relation between friction factor (f) and Reynolds number 
(Re) which follows 0 

P 4 = ~ + 0.55 fC 5 (2.8) 

which is similar to eqn. (2.6) obtained by Ahmed & Sundada 
previsouly, except for 0.55V"0^ m place of 1.0 where C ~ is a 
variable depending on size, shape, surface roughness, porosity. 

Todd and Tyagi (1970) obtained values of a and b of 
Porchheimers equation and based on the data a graph between 
friction factor and Reynolds number was presented which 
Cont-radicts the findings of Ahmed and Sundada (1969) » Werd 
(1970), Rang arm cl ha Rao and Suresh(1970). 

Izbash and leleeve (1971) has defined two limits 
of Reynolds number R e q im first boundary Reynolds number, 
above which results obtained by assuming pure laminar flow 
conditions leads to an error of 15 per cent. 

Pe 4rrfl - — s&eond Imamdary— tReyn^idrS-- nussber-, — abave—whictu re-su-lts 
obtained- by using square law leads to e rror less th an 45 -per~ceu%-. 
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II 

Re lim secon< 3- boundary Reynolds number, above which results 
obtained by using square law leads to error less than 15 
per cent. 

With these definitions new laws for heads loss have 
been defined 


I 


I 


I 

lim 


II 

lim 


= 216 


= 216 


Re f_ (J=S)5 


gd" 


'lim v n 


C d 
o 


fT) II \2 ,1-n \3 
J (Re lim } ( “1T ) 


(2.9) 

( 2 . 10 ) 


where 


I = hydraulic gradient 

L< = dynamic viscosity 

g = acceleration due to gravity 

n = porosity 

d = size of gravel 

C 0 = 1/s, where s is coefficient of microroughness 


P = coeCPicient depends on cross-sectional shape. 


Subramanya and Rao (1971) obtained relationship, 
based on experimental work between Rorchheimer 's coefficients 
a and b and the specific per-meability (hs) which are - 


a Ks 



( 2 . 11 ) 


= 1.42 x 10" 4 


b Ks 


( 2 . 12 ) 
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The specific permeability (ks) was obtained by 
ks = (k/y ) (2.13) 

\jtsd.0S'iy 

where y is d ensity of water and 
k is the permeability. 

The above two relations are valid m the range of 
Reynolds number from 10 to 200 which is one of the limitations 
to use the relations. Since size shape, surface roughness 
and porosity, which are important factors in the flow of 
water through porous media are not included m the relation- 
ship discussed, it was felt necessary to study these parameters 
and to get a relationship between them. 



CHAPTER 3 


STUDY O F ITONLIUCAR FLOW PARAPUTERS 

3.1 &ENERAL 

(U tbvxtU^' 

It is clear from the review in chapter 2, that even- 
though some studies hove been made on the non Darcy flow 

L YA 

through soil, no effort seeras to have been devoted w relating 
the coefficient a and b, especially the later, of the 
Eorchheimer 's equation, with the physical characteristics of 
the porous media . The available data is insufficient m 
detail to establish a relation between the parameters of 
interest and hence, this study. 

3 . 2 E XPERIMENTAL SS T-UP ; 

A detailed diagram of the experimental setup is shown 
m Dig. 3.1. Continuous supply of water is provided by the 
overhead tank which maintains a constant head, the permeameber 
is 15 cm m dia. and 37.5 cm long with a number of pressure 
tappings on it. These tappings are located m tw6 rows 90° 
apart along the length of the permeameter and are staggered. 
This arrangement of pressure tappings gives a continuous 
measurement of the hydraulic gradient along the length of 
the permeameter. Any local variations or inadmissible values 
are easily located. A perforated plate is placed at the 
bottom of the permeameter to support the material. A valve 
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is provided at the junction of the permeameter and the delivery 
pipe to regulate the flour. The delivery pipe, which has a 
slit m the base rests m a measuring tank of size 90 cm x 75 cm. 
This measuring tank is fitted with a 90° v-notch to measure 
the discharge. A dram adjacent to the measuring tank 
collects the w^ter which m turn flows to a sump and a centri- 
fugal pump lifts the water to the overhead tank to complete 
the cycle. The measuring tank is provided with a device to 
read the water-level m it. Pressure difference between any 
two tappings is read on a mercury manometer. A wooden mesh 
is provided near the delivery pipe m the measuring tank to 
dissipate all the energy of che discharged water. 

Material for investigation has been selected from 
-the gra /el > and glass balls. The sizes of gravel are 
0.318, 0,636, 1.115, 1.75, 2.38, 3.33 and 4.62 cms. The 

sizes of glass balls are 1.3, 1.7, 2.0 and 2.5 cms. Material 
used is uniform m size. 

It was decided to test each size of the material at 
three different porosities, very dense, very loose and a 
porosity m between the two extremes. 

3 . 3 PACKI NG OP THE MAT GRI Aly 

Small sizes oC the material were packed by a hapmer 
(5 kg.) to get the required porosity. The bigger sizes (_> 1.5 cm) 
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were packed by hand to get the required porosity. While 
packing the material , following precautions were taken. 

(1) Material was packed uniformly 

(2) Top of the packing was made level and average 
depth was measured to get the height of packing. 

(3) If hamper has been used to compact the material 
than care was taken to see that the material w as 
not crushed by the hammer blows. 

(4) If material has been packed loosly by hand than 
proper interlocking of individual gravel was 
assured. 

(5) All attempts have been made to get three distinct 
por^sitie s. 

3 . 4 PREPARATI ON 02 IH 'J St lTUP; 

The permeameters filled with the material to be 
tested was fixed m position. Both top and bottom joints 
were made leak proof:. TT ercury levels m the manometer were 
checked initially and scales adjusted to read similar 
readings m both the limbs. All the pressure tappings 
were closed except the one which has been connected to the 
nanometer. The by pass valve of the manometer was opened. 
The initial reading of waler level m the measuring tank 
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was noted. The pump was started and the water was allowed 
to flow through the permeameter. In the beginning it 
was noticed that air bubbles were coming out from the tube. 

As the bubbling stopped the other tube was connected to the 
required pressure tapp mg . All the connections were checked 
and assured to be leakproof. 

3 . 5 METHOD OF OBSERV ATIONS • 

Observation were taken for decreasing and increasing 
flow rates. For any flow rate 9 manometer readings we re 
taken for 9 different combinations of the pressure tappings. 
The water level in the measuring tank was noted. Thus 9 
pressure tapping readings and one water level reading 
completes one set of observations for one flow rate. With the 
help of regulating valve flow was slightly reduced such 
that water level m Lhc measuring tank was reduced by 
approximately 2 eras. I non another set of observation was 
recorded. Thus 5 sets were taken for decreasing flow rates 
and 4 sets for increasing flow rates. At the end, the 
tenfperature of the water was recorded. Each test was repeated 
till reproducible results were obtained. 

During observations the following precautions were 

taken: 

1 . If mercury level m the manometer was found to be 

oscillating, the average value was recorded. 
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2, By-pass valve of the manorneter was operated very 

carefully* It is necessary to do so, specially m 
the case of lo v v pressure differences. 

3* All efforts were made to dissipate all the energy 

of water m the measuring tank. This helps m 
reading accurately the water level. 

4. Bottom supporting inesh of the permeameter was 
replaced by a coarser mesh before testing larger 
sites of the gravel. 

5. Only silt free clear water was allowed to flow through 
the permeameter. 

3 . 6 RES ULT S AND DISCUS S ION 

Prom the observations taken, the discharge 
m cc/sec., velocity (v) l o cn/sec . and the average 
gradient (l) were obtained. Then the values of i/v, 

Reynold's number (Re) and friction factor (P) were also 
cal culated. 

Coefficient a end b were obtained by a graphical 
method and were verified by a statistical fit. 

Graphical Method 

Points were plotted on a graph with i/v and v taken 
on vertical and horizontal axes respectively. These points 
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were found to follow a straight line. The intercept of 
this straight line on tie vertical axes and its slope 
(multiplied by scale ratio) give the coefficients 
a and b respectively. 

Statistical Fit 

A. computer programme was prepared to fit a straight 
line to the observed data using the method of least squares , 
and the two parameters a and b obtained. 

Experimental v clues of 'a' and 'b' as obtained are 
shown m Table (3.1). 

The variation of l/v and v is shown m Fig. 3.4 
through 3.10 for the different sizes of gravel and Fig. 3.11 
for all sizes of glass balls. These plots confirm the vali 
dity of the Forcbheiraer *s equation for the ranges of 
hydraulic gradients and velocities used. The flow can be 
classified as 'Steady Inertial' as per Wright (1968). 

The linear fit shown is extended to meet the Y-axis but 
it is believed that the actual curve would be different 
from this extended line due to the existence of laminar flow 
regime where Forcv’s law is valid. 

The coefficient 'a' obtained by extending the fit 
straight line fit would bo slightly smaller than the actual 
value. It is obvious that uhe values of 'a' for each sice, 
increases with decreasing porosity 'n' (or void ratio-e). 



24 


This is a consistent result as 'a' is equivalent to 
the inverse of the coefficient of permeability and the 
lat er should decrease ' r ith decrease m 'n' or ’e 1 . 
Jagadeesha Rao (1971) attempted to corretate r a' with the 
media characteristics and gave the following” 


’a * 


= 10 ' 


-5 


(3*1 ) 


where K is the specific permeability, 
s 

The nonlinear parameter 'b ' increases with 
decreasing values of porosity for all sizes of gravel and 
glass balls. For the smallest sue of gravel used (size 
0.318 cm), the porosity could not be varied much and hence 
the slope of i/v vs v curves varied by small amount only. 

For sizes which could be tested over larger ranges of porosity 
■foe slope and hence the coefficient ! b' varied over a 
greater range. 


By analogy with pipe flow, the inertial loss can 
be written a s 


2 

f L 1 v" 

2g D-j 



f v 
2g D 


(3.2) 


where h^ is he- 1 ! loss, f is a constant depending on 
Reynolds number, and are the length and dia of 

flow and v is the velocity. Comparing this with second 
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* 2 

term of the Porchheimcr * s equation (i = av + bv ), one 
can write 


b = 


cons bent 



(3.3) 


where D-| is the size of capillary of the porous media. 
Since this 3^ decreases with decreasing porosity and 
particle size, one could anticipate to g et higher values 
of ' b'. The same feature is observed m all the c?ses. 

Pig. 3.11 depicts the results of all sizes of glass 
balls tested. With the gloss balls it was not possible to 
get more than two porosities for each size. However, the 
trends are very smaller to the observed with gravel. 

The dotted lines show the Reynolds number contours 
on some of the curves drawn for the sizes 0.318, 1.75 and 
4.62 cm of the gravel. These contou.s were drawn by first 
writing the corresponding values of the Reynolds number 
on each point then lniernolvating. It has been found that 
the contours are nearly vortical m the case of 0.318 cm 
size and inclined towards right side with decreasing 
slope for increasing discharge and the size of the gravel. 


The influence or bp a pore channel size which itself 
depends on particle size, shape and porosity of the medium, 
onthe nonlinear parameter is studied further. If the 
variations m ’b’ are related to the particle size, the 
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observation made earlier that ’ b * decreases with an 
increasing particle size and porosity, gets confirmed. 

The available data (Coeff. b) from Dudgeon ( 1964 ), 
Ahmed and Sundada (1969), Ranganadha Rao and Suresh (1970) 
and the data of this work have been plotted with the size 
of the gravel on a log-log plot (Dig. 3.12). It has been 
found that all points lie closely along a straight line. 

The equation of the straight line was found to be 

log^ = 0.712 Jjog lQ - 1.237 (3.4) 

equation (3.4) is further simplified to 

b = -§5 (3.5) 

d 1 

The equation 3. / ! or 3.5 does not include the 
parameters other than size. Hence this equation gives 
only the approximate variation of coefficient 'b ' with the 
size of the grovel. However, Fig. 3.13 shows the same vari- 
ation on natural scale. 

Fig. 3.14 shows the variation of ’b' with porosity 
for different sizes of gravel and glass balls. The effect 
of increasing b-value with decreasing porosity is obvious. 

An important aspect 1 . 0 . the shape and the roughness of 
the pore channel can be seen from this figure. In case 
of glass balls the channels are bound by smooth spherical 
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surfaces while mease ol gravel the particles are 
angular ^nd their surface is rough. Between the two factors 
shape and surface roughness, the later is supposed to have 
a dominating influence. The 'b ! values are smaller for 
glass balls as their surfaces are relatively smooth. 

Coefficient 'b ! has been related to the parameters 
like - equivalent sphei ical dia 'ds', porosity 'n' , 
surface roughness and shape of gravel and glass balls. An 
empencal formula relatin' the above parameters and the 
coefficient f b f was fount) to be 


where R’ is a parameter for the steady inertial flows. 
Since it was not possible to measure surface roughness 
and shape factor separately their effects have been 
included m one parameter R’. The experimental values of 
R 1 for angular and highly rough surfaces like gravel is 1.0 
and for spherical and smooth surfaces like glass balls is 
1.25. Therefore for intermediate rough surfaces and shapes 
like river gravels, the parameter R’ may be interpolated. 

A similar emponcal formula has been suggested by Aravin 
(1965), m which power of (ds n) is 2.0 for turbulent 
flows, irrespective of shape and roughness of the particles. 
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A new non-dimensional parameter has been obtained 
by multiplying the coe evident f b ' by the acceleration due 
to gravity (g) and the equivalent spherical dia (ds). 

This new parameter has been studied m relation with the other 
non-dimensional parameter s i.e. porosity ’n’. Pig. 3.15 
shows the variation oC b.g.ds with'n which also clearly 
dipicts the influence of surface roughness and shape of 
gravel on value of ’b’. Hence this graph can be used to 

predict ! b* for the hnoon values of porosity and the 
equivalent spherical dia (ds). 

Pi . 3.16 illustrates the P vs Re variation. 

In this graph the points are found to be spread m a small 
range. This spread is due to the large variation m the 
porosities m which the each material has been tested graph 
shows two curves i.e. one for gravel and another for glass 
balls. However 5 the trend of the graph is similar to the 
graph obtained by the other investigator. Curve is nonlinc >r 
for Re greater than 8.0. 
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TA3L C 3.1 



SUPER II ISill AS VALUE'S OP ’a.’ and__ 

•b ’ 





Type of 
Material 

Nominal 

size 

(cm) 

percentage 

porosity 

’a’ 

Sec/cm. 

•b 1 

Sec 2 /cm 2 

Gravel 

0.318 

42.0 

0.288 

0.093 



41.6 

0.300 

0.101 



41.2 

0.300 

0.103 



40.2 

0.310 

0.107 

t ? 

0.638 

43.5 

0.06 

0.042 



40.8 

0.08 

0.048 



36.0 

0.10 

0.067 

T ? 

1.115 

43.0 

0.016 

0.026 



38.0 

0.076 

0.041 



34.0 

0. 160 

0.054 

t f 

1.750 

46.5 

0.01 

0.0102 



41.5 

0.02 

0.0105 



36 . 0 

0.035 

0.0173 

f f 

2.38 

44.7 

0.005 

0.004 



40.5 

0.005 

0.00824 



35.5 

0.007 

0.0148 

t t 

3.330 

50,0 

0.008 

0.0021 



46.6 

0.01 

0.0029 



43.0 

0.055 

0.0040 




Gontd. next page 
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Type of 

Nominal 

percentage 

V 

’b ' 

Material 

size 

(cm) 

Porosity 

Sec/cm 

Sec 2 /cm 2 

G-ravel 

4.62 

50.0 

0.002 

0.0010 



46.5 

0.004 

0.0019 



41.6 

0.028 

0.00372 

Glass Balls 

1.30 

35.5 

0.0324 

0.0081 

t ? 

1.70 

36.7 

0.0133 

0.0066 



36.1 

0.0293 

0.00677 

1 f 

2.00 

38.3 

0.0192 

0.00448 



36.7 

0.0278 

0.00525 

t t 

2.50 

39.2 

0.0140 

0.00275 



38.1 

0.0209 

0.00347 



Over head 
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CHAPTQR 4 


D EVELOPMENT OP A C0KTIC4.L P JRKEMSTER 

4. 1 GENERAL 

In many field situations, flow may be convergent 
with large variation in Hydraulic gradient over a short 
distance, e.g. flow in to a well, rockfill and earth dams. 

It has been established tb ,t most coarse grained soils exhibx 
nonlinear flow characteristics. The flow parameter ’K r , 
the coefficient of permeability will not be a constant 
over the total range of hydraulic gradients. In establishing 
discharges for such situations, often, permeability is det- 
ermined using a setup bl>- t essentially devel^s a constant 
or nearly constant hydraulic gr-dient. 

Therefore, it is fait necessary to develop a 
new permeameter that con generate variable hydraulic 
gradients and to test the soils to get some average 
coefficients from the s^ie. A conical permeameter is 
constructed, necessary formula for flow through the same 
are derived and some preliminary tests were carried out. 
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4.2 DERIVATION OF THE FORT TJLAE : 

Formulae have been derived for the discharge 
through the conical permeameter using: 

(a) Darcy's la?/ and (b) Forchheimer ' s equation. 


4.2.1 Derivation Using Darcy's law:- 

OASB 1 : Constant head case (Fig. 4.1) 

r.j and rg are the bottom and the top radii of the 
conical permeameter . 

1 is the height of permeameter 

s is the taper m the sides of the permeameter 

h is the head difference between inlet and 

exist at any instant. 

K is the coefficient of permeability, 

q is the discharge. 

Consider an element dx 3 at a distance x from bottom; 
m the permeameter. From Darcy's law 


K i x A x 


or 


3L 


KA. 


x 


(4.1) 


where i and A „ are the nydraulic gradient and cross- 

X X. 

section area of the permeameter respectively at the strip 


m consideration 
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2 2 
= ic r = n; (r. + sx) 

•A. X I 


X 


ah 

dx 


Km(r^+ s x) 1 


(4.2) 


or an = a- 

Kte 


dx 

(r, + s x)' 


(4.3) 


Integrating equation (4.3) and taking limits 


h = 


ll 

Km r 2 


(4.4) 


or 


1 = 


h K % r 1 r 2 


(4.5) 


CASE 2 : Variable heed case. 

II A is the cross-sectional area of the inlet 
tube, then 

Equations (4.5) and (4.6) 

h r i r 2 IC ” _ , an (4.7) 

— - - A dt 

Let h change from to E 2 when time changes from t^ 
to 1 2 * 
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(4.7) 


Separating the variables and integrating equation 


r 0 r, Ktc ^2 

_2_J — j at 

la t. 


K 

f 


1 


r 


2 


dh 

h 


or 


T i ^ 

I A 


(t 2 - t 1 ) 


In 



) 


(4.8) 


(4.9) 


Equation (4.9) gives the required relationship. Howeger, 
if r = r^i = r 2 is substicuted m it. 


rr 2 II . 

IJLi™ = i n ( -1-) (4.10) 

1 1 X 2 

which is the formula for cylindrical permeameter. 


4.2.2 : Derivation Using Forchheimer ' s Equation: 

CASE 3 .: Constant Head Case (Eig. 4.1). 

dotation remains some. Considering the element' 
dx at distance x Cron bottom. Applying Forchheimer T s 
equation to bins client. 


'x 


= av,, + b v. 


(where v = 

X. 



(4.11) 


or 


A + b A 


a'T 


( 4 . 12 ) 


4 
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7t(r^ + sx)' 


Tc 2 * (r^ + sx) 4 


( 4 . 15 ) 


Separating the variables and talcing the limits of 
equation (4.13). 


h 

/ dh 
0 


f a . 1 . dx _ 

0 Tc(r^ + sx) 1 


L . 2 , 

f b q dx 

0 ^ 2 (r 1 + sx) 4 


( 4 . 14 ) 


Integrating and simplifying equation ( 4 . 14 ) 


a q 1 


% r^ r 2 


, 2 X / 2 2 \ 
o q L (r 1 + r ? + r 1 r«) 


J v.- 1 - -j ■ 

3 n: 2 r; 


( 4 . 15 ) 


CASE 4 : The Variable Head Case. 


Rewriting equation (4.15) as 


K T , 2 ^ 2 n 2 T 22 

b L (r^ + r 2 + r 1 r 2 ) q + 3 alixr^r 2 q 


- r^ r| h = 0 


( 4 . 16 ) 


2 2 

and substituting P = b 1 (r^ + r 2 + r.j r 2 ) 

2 2 

n = 3 q L 71 r^ r 2 

2 3 3 

X = 3 x rj r 2 


Milt 
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Equation ( 4 . 16) reduces to 


Pq +Qq-bX = 0 


(4.17) 


Roots of equation (4.17) are 


q = 


Q + t T + 4 P h X 


2P 


( 4 . 18 ) 


Since q can not be negative. 


q = T + 4 P X h 

2 P 


(4.19) 


In variable head case 


« dh 
q A dt 


( 4 . 20 ) 


Equating (4.19) and (4.20) 


Adh _ 
dt 


Q + { Q 2 + 4 P X h 
2 P 


( 4 . 21 ) 


let h change from H„ to H 0 when the chance m time is 


L 1 


from t.j to tg. 


Separating the variables and integrating Eq. 4.21 


H. 


I 

Er 


A dh 

-Q + f q 2 + 4 P X h 


f 

t 


dt 

2P 


( 4 . 22 ) 
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On integrating and simplifying equation (4.22) 


(i/X) V Q 2 + 4- P X H. - (JL/X) Y Q 2 + 4 P X E, 


+ (AQ/X) In rf Q 2 - 4 P X H, - Q) 


In ( Y Q 2 + 4 P X H 2 - Q) = tg - t 1 


(4.23) 


Equation (4.23) further simplifies to 


In 


4 P X H 


1 + 


1 


- 1 


Q‘ 


I - \ ± JL X H 2 " - 1 


l 1 + p- 

' Q 2 


J 



4 P X H 


(4.24) 


1 


) of equation (4.24) on 


The expression (1 + -p- 

Q 

substituting bach the values of P ? Q and X and simplifying. 


1 + 4 P X H, 


Q £ 


A. u 

2 H 1 


a 


1 + B 


(4.25) 
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where . r. r Q 

t, i~ ( _1 + -£ + i ) 

B ~ 31 ^ r 2 ^ u 


(4. 26) 


The equation (4.24) reduces to 


B b H 


1 + 


1 


Z 


lu 


a 


>/ 


/ B b H, 


1 + 


a 


7i r 1 r 2 
11a 


- (/ 1 + 
I 


B b H 1 ; B b H 2 

SL 


(4.27) 


If (b/a^) is much less than unity, expanding 
B b EL 1/2 

(1 + - 2 — ) bmorialy a nd considering the first two 


a 


terns only 


( 1 + 


b H 1 1/2 = 1 + B b H 1 


a 


"2““ 


) 


2a 


^ + Higher power 


terms 


(4.28) 


Substituting the quantities of (4.28) m equation 
(4.27) and simplifying. 


In 


K, 


% r^ r 2 At 

1 A a 


B b 
2 a 2 


(H 1 - H 2 ) 


(4.29) 
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Similar equation can be formed for 2nd and 3rd 
values of Hg and H^» with 'die condition that Alt remains 
same . 


H 0 tt r, r 0 a t 

In ( J- ) = 1 2 


B b 


Eg 


LA a 


2 a 


2 (h 2 - h 3 ) 


(4.30) 


Subtracting (4.30) from (4.29) and simplifying 


In ( 


H 


H 1 H 3 


■) = 


2 of 


( H 1 + H - 2H 2 ) (4.31) 


Since there are two unkno wns ' a' and ‘ b the two 
equation (4.29) and (4.31) can be used to determine them. 

If the value of nonlinear coefficient T b’ is set 
zero in the eq_uabion (4.29) , the equation reduces to the 
equation (4.9) derived from the Darcy's law. 


4 . 3 EXPE RIMENTAL - 70RK : 

Coefficient of permeability (If) has been obtained 
to investigate the effect of variable gradient on it with 
the help of equation (4.9). To do this three values of 
’K' have been obtained, corresponding to three porosities 
for each type of permeemeter and for sand. 
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The setup is shown m Pig. 4.1. It consist of 
a permeameter (conical or cylindrical) and a glass jar. 

Por filling the permeameter, required volume and hence the 
required weight of clean and dry sand was obtained for a 
particular porosity. This sand was poured- m such that it 
completely fills the permeameter uniformly. After this 
cover w^s bolted assuring no leakage. The inlet tube of 
the permeameter was then connected to the jar. 

Clear water was now poured m to the jar up to the 
level slightly higher (2-5 eras) than the mark at which the 
observation was to be recorded. Out let valve was opened 
and as soon the water level m the jar crosses the pre- 
decided mark, the stop watch was started and the heights 
were recorded after eveiy minute. 

Before taking the observations, the sand was 
saturated by allowing the water to flow through it for at 
least 15 minutes. Whole experiment was repeated untill 
reporducible results were obtained. 

4.4 RESULTS 

Prom the recorded observations, the values of 'K ' 
were obtained by using equation ( 4 . 9 ). Thus the values 
of K, obtained from cylindrical and conical permeameters 
with the contraction ra'cio of 1:2 and 1:4, were plotted m 
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a single plot (Fig. 4.2). A single curve fits all the 
points 5 thus establish ang that for the range of hydraulic 
gradients existing m the permeameter, K ti* coefficient of 
permeability is nearly independent of the gradient. 

Some tests have been conducted and analysis made usin' 
equations (4.29) and (4.31). The results are in conclusive 
and hence are not presented here. 




DIAMETER RA 1 



VOID RATIO 



CKAPTCR 5 


cQTTCLUDr'r Rrmpxs 

In this investigation, the flow through soil 
(porous media) is studied m the 'Steady Inertial Regime’ 
wherein deviations from Darcy’s law are ooserved, but 
the flow is still lamn^r. To characterise the velocity 
hydraulic gradient relscion, the Dorchhe inter equation 
( i = av + bv ) is ch >osen m preference to the Missbach 
equation (l = cv 1 ) as the former is derivable (reducible) 
from the Navier - Stoke 1 s eouations. 

The variations m I'he values of the coefficient 
' b f have been studied lfch sice and porosity individually. 
The effect of surface roughness end shape of the grains 
could not be seperated -nd hence are studied in a combined 
form. It is found that ’b* decreases with increasing 
size as well as porosity and the variations between them 
have been shown m fig. j>.12 and 3.14. These graphs also 
depict the effect of surface roughness and shape of gravel. 
A rel^rion-chip has b^en developed between coefficient ’b' 
and the sice (d) of gravel, which follows: 


0.0185 



b 


(5.1) 
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From the above equation an approximate estimate 
of the 'b' value may be obtained for a given size. 

Another relationship has been proposed for calculating 
a fairly accurate value of the coefficient ! b' m which 
the effect of size (ds), porosity (n), shape and surface 
roughness (R) of tne gravel have been combined together. 
The relationship is as follows 

b = — ( 5> 2) 

(ds n) R' 

In the above equation parameter R' is 1.0 for 
angular and rough surfaces end is 1.25 for spherical 
and smooth surfaces. A new nondimensional parameter 
(h.g.ds) has been obtained and its variation with 
porosity has been shown m Fig. 3.15. This graph also 
shows the effect of shape and roughness of the gravel 
and may be used to ascertain the coefficient * b f . 
Eventhough a relation apeears to exist between the 
non-dimensional parameter (b.g.ds) with porosity (n), the 
points obtained are few to arrive at Che relationship. 

A large number of tests aced to be carried out further. 

A curve (Fig. 3. IS'' between friction factor and 
Reynold’s number has be >n drawn m which points are 
found to be scattered. The reason for such scattering 
of points is the range of porosity m which each size has 
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been tested. This is also obvious from the fact that for 
different materials at a parcicular Reynold's number (Re), 
friction factor (IF) can not be the same since it is a 
function of particle sire, shape, surface roughness and 
porosity. 

In any field situation, the hydraulic gradient is not 
constant but varies in sp~ce the coarse grained soils 
(like sands, gravels etc.) exhibit nonlinear relation 
between hydraulic gradient and velocity. So it was felt 
necessary to develop a permeameter which can generate a 
range of hydraulic gradients and get an average coefficient 
of permeability for th? range generated conical permeameter 
is one such which satisfies the above requ trement . 
Expressions have been derived for head loss and discharge 
through the conical penaeaneter for either linear 
(Darcy’s equation) or non linear (Eorchheimer ’s equation) 
flow conditions. Some typical results for K are obtained 
for fine sand which compare with those obtained using 
a cylindrical permeameter. 
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